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Abstract 

In the scenario where the dark matter (DM) particles XX annihilate through 
a resonance particle R, the constraint from DM relic density makes the correspond¬ 
ing cross section for DM-nuclei elastic scattering extremely small, and can be below 
the neutrino background induced by the coherent neutrino-nuclei scattering, which 
makes the DM particle beyond the reach of the conventional DM direct detection 
experiments. We present an improved analytical calculation of the DM relic den¬ 
sity in the case of resonant DM annihilation for s- and p-wave cases and invesitgate 
the condition for the DM-nuclei scattering cross section to be above the neutrino 
background. We show that in Higgs-portal type models, for DM particles with 
s-wave annihilation, the spin-independent DM-nucleus scattering cross section is 
proportional to the ratio of the decay width and the mass of R. For a 

typical DM particle mass 50 GeV, the condition leads to TR/mR > 0(10 ^). In 
p-wave annihilation case, the spin-independent scattering cross section is insensi¬ 
tive to VR/niR, and is always above the neutrino background, as long as the DM 
particle is lighter than the top quark. The real singlet DM model is discussed as a 
concrete example. 
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1 Introduction 


Dark matter (DM) contributes to 26.8% of the total energy density of the Uni¬ 
verse [1], yet its particle nature remains largely unknown. The leading candidates for 
DM are weakly interacting massive particles (WIMPs). WIMPs can naturally obtain the 
observed relic density, and the predicted cross sections of the WIMP-nuclei scattering 
are usually within the reach of the current DM direct detection experiments. In the case 
where the DM annihilation cross section times the relative velocity crUrei is a constant, 
such as that in the simple s-wave annihilation cases, the DM relic density can be calcu¬ 
lated analytically using the standard freeze-out approximation. The connection between 
the DM relic density and the DM-nuclei scattering cross section can be straightforwardly 
established. 

However, in many DM models and DM interaction mechanisms, the velocity depen¬ 
dence of crUrei can be complicated. For instance, to explain both the relic density and 
the cosmic-ray positron excess observed by PAMELE [2], Fermi-LAT [3], and AMS- 
02 [4,5], the mechanism of Sommerfeld enhancement has been invoked, which introduces 
a velocity-dependent DM annihilation cross section [6-13] to account for the larger cross 
section required by the data [14-16]. 

In a wide class of DM models, the DM particle y can annihilate into the standard 
model (SM) particles through an s-channel resonance particle R. Such as the singlet 
scalar DM models [17-26] , the left-right symmetric models with extended stable scalar 
sectors [27-34] and the fermionic DM models [35-41]. The presence of the scalar R can 
also play an important role in electroweak phase transion [24,39,42-50] and modify the 
interpretation of the DM-nuclei scattering [51,52]. 

Near the resonance point the kinetic energy of the DM particles is non-negligible, 
which makes aUrei velocity dependent, and leads to the enhancement of DM annihilation 
cross section at lower tempertures, the so called Breit-Wigner enhancement [53,54]. In 
the scenario of resonant dark matter (DM) annihilation, under the constraint of DM 
relic density, the cross section for DM-nuclei elastic scattering can be extremely small 
such that it can fall below the background induced by the coherent neutrino-nuclei 
scattering, which make it undetectable by the current DM direct detection technology. 
It is of importance to know under what condition this phenomena will occur. In order 
to establish the correlation between the DM relic density and the DM-nuclei scattering 
cross section, it is useful to have analytical expressions for both quantities, which is 
however challenging, due to the complicated velocity dependence of the DM annihilation 
cross section in the case with resonance. 

The Boltzmann equation which governs the evolution of the DM number density is 
usually solved analytically by using the standard freeze-out approximation [55]. However, 
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when DM annihilation takes place near a pole in the cross section, we cannot use the 
standard method as (crUrei) does not have a simple analytical form [56]. In Ref. [57], it 
was proposed to analytically calculate (crUrei) using the (5-function approximation, if the 
resonance has a very narrow decay width. But this method fails in the case where the 
DM mass is greater than a half of the mass of the resonant particle, namely, above the 
resonance. 

In this work, we present an improved analytical calculation of the DM relic density 
in the case of resonant DM annihilation for s- and p-wave cases and investigate the 
condition for the DM-nuclei scattering cross section to be above the neutrino background. 
We show that in Higgs-portal type models, for DM particles with s-wave annihilation, 
the spin-independent DM-nucleus scattering cross section is proportional to Tr/itiji, 
the ratio of the decay width and the mass of R. For a typical DM particle mass ~ 
50 GeV, the condition leads to Tr/itir > 0(10“^). In p-wave annihilation case, the 
spin-independent scattering cross section is insensitive to TR/mR, and is always above 
the neutrino background, as long as the DM particle is lighter than the top quark. 
As an example, we calculate the spin-independent cross section both analytically and 
numerically in the real singlet DM model with resonant annihilation. We show that the 
predicted cross section in this model is always above the neutrino background. In order 
to cover the full parameter space of this model, the required sensitivity should reach 
~ 2.6 X cm^ for the next generation direct detection experiments. 

This paper is organized as follows: In Sec. 2, we outline the freeze-out approximation 
of Boltzmann equation, and propose an approximate formula for the relic density of a 
scalar or fermion DM particle which annihilates through an s-channel scalar resonance 
which has a narrow decay width. In Sec. 3, we analyse the constraint of neutrino back- 
groud has on DM direct detection experiments on the resonance point. In Sec. 4, we 
analyse the direct detection of the real singlet DM. Some discussions and conclusions are 
given in Sec. 5. 

2 DM relic density from resonant annihilation 

The time evolution of the number density n of the DM particle y is described by the 
Boltzmann equation 

^ = -3Hn - {avrei)in‘^ - nl^), (1) 

where Ueq is the equilibrium number density of y, H is the Hubble parameter, and ((JUrei) 
is the thermal average of the total annihilation cross section times the relative velocity 
frei of the annihilating particles. In the non-relativistic case, the thermally averaged 
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cross section can be written as 

^ 3/2 roo 

(aVrel) = J ( 2 ) 

where x = rn^/T with T the temperature of the photon in equlibrium and the mass 
of the DM particle. Dehning Y = n/s as the comoving density of particle x with s the 
entropy density, Eq. (1) can be rewritten as 


^ _ I Tig* Mpim^{avrei) ,2 

dx V 45 

where Mpi = 1.2211 x 10^® GeV is the Plank mass, and 
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where g^s and hea are the effective relativistic degrees of freedom for entropy and energy 
density, and 
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where g is the internal degrees of freedom of the DM particle y. The decoupling tem¬ 
perature X/ is dehned as the temperature at which the DM particles start to depart 
from the thermal equilibrium, and the density Y is related to the equilibrium density 
Teq by Y{xf) = (1 -|- c)Yeq{xf), where c is a constant of order unity. The value of x/ is 
approximately given by [58] 


Xf ^ \n[0.038c{c + 2)Mpim^ggJ^^(av.ei)] 

~ lnln[0.038c(c 2)Mpim^gg~g^^{aVrei)]- (6) 

The value of c is usually taken to be one, which leads to a good £t to the numerical 
solutions of the Boltzmann equation. The DM number density in the present day Yq can 
be obtained by integrating Eq. (3) with respect to x in the region Xf < x < 00 , 
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the function J/ is dehned as 
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where we have used the definition of (aurei) in Eq- (2). Exchanging the order of the 
integration in Eq. ( 8 ), Jf can be represented as [56] 


Jf = 
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rel" 
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(anrei)nreierfc( y^ni.el/2)dnrel- 


The relic density of y is obtained from Yq as 
= 2.755 X lO^Er 




45 1 

-Kg^ MpiJf 


-GeV 
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2.1 The case of s-wave annihilation 


We first consider a real scalar DM particle x annihilating into SM particles through 
exchanging a mediator particle R in s-channel. The interaction between x and R can be 
written as £ D ^nxxR + where /i is a dimensional coupling constant. The 

term fiuRH^H leads to a mixing between R and the SM Higgs boson H. x can be 
stable due to a Z 2 symmetry x XX —x- The annihilation proceeds through s-wave, the 
corresponding cross section multiplied by Urei is given by 

2 ^^ £. r(fi- ^ V,) 

(s - m |)2 + m|r| 2 m^ 

where and T/j are the mass and total decay width of the resonance i?, s is the 
Mandelstam variable, in the non-relativistic case s ~ 4m^ + R* -x Xi stands for 

any possible decay mode of R* and '^^^{R* -x X^) is its total decay width. If is 

close to the resonant point (-y/s ~ 2 m^), ^{R* —t X^) can be taken as the total 

decay width T^^ and Eq. (11) can be rewritten as 
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From Eq. ( 8 ) and Eq. (12), the expression of can be rewritten as 

7H'^^rel 




Jf= / ^erfc(vT7nrei/2)-^^—^ 
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There is no analytical expression available for Jf. If -C 1 and -C using 

the relation 


lim 


7h 


7 h^o (n2^j - e/j)2 + y| 

the value of Jj can be approximated as [57] 


= - ^r), 


Jf 


TTI 

Jj = j^erfc(^/x/€H/2), for £r > 0 and 7 r <£ 1. 


(15) 


(16) 


Note, however that this approximation is only valid for en > 0. 

In this paper we present an improved method to evaluate Jf which is valid for both 
Cij > 0 and e/j < 0 with a reasonable precision. If > 0 and y^j 1, the integral 
of Eq. (14) dominates in the narrow region near the point Urei = \/^- In this re¬ 
gion the complementary error function erfc(y7i/nrei/2) changes very little. We can take 
erfc(y^frei/2) ~ eTfc{^Xfeji/2), therefore 


Jf 
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—h arctan — , for ej? > 0. 

2 IR, 


(17) 


For eR 3> yij, arctan(eij/yi{) 7r/2 and Jj ~ J^. Likewise, if e^j < 0 and the absolute 
value of eR approaches zero, the integral of Eq. (14) dominates in the region near Urei = 0 
and we can take erfc(y^r’rei/2) ~ 1, then 

,2 
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ja _ h' 

^ Amt 
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- + arctan^|, for £„ < 0. 
2 7h, 


(18) 


In Figure 1(a), we show the differences in the approximate analytical results Jf, Jf, and 
the numerical result JJ in a specihc case where the parameters are taken as /i = 1 GeV, 
m_R=200 GeV, = 10, Xf = 20, and Tr = 0.001 GeV. From the figure, if < 
100 GeV, the analytical result JJ agrees with the numerical result JJ very well, the 
relative error is less than 2%, and the analytical result Jj can obtain the same precision 
if cr S> 'Jr is satished. If > lOOGeV, the approximation of Jj is no longer valid, but 
JJ still agrees with the numerical result well near the resonance point with the relative 
error is within 12% (in the region 100 GeV < < 100.01 GeV). The error increases 

with leaving away from the resonance point. 

From Eq. (10), (17) and (18), the relic density of y can be represented as 
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for eR > 0. 
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for Sr < 0. (19) 
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On the resonance point {m^ ~ 771^/2), we find 


n o / 45 8m‘l ^ 

nh^ ^ 2.755 X 10%/ - ^^GeV\ (20) 

V Tig, vr/iWpi ^ ^ 

Eq. (20) shows that the relic density is not sensitive to the decay width of R on the 
resonance point. 

If X is a complex scalar, the interaction between y and i? is £ D fi\x\‘^R, the expression 
for crWrei and relic density is identical to the case of real scalar DM. 


2.2 The case of p-wave annihilation 


If X is a Dirac particle, the interaction between x and R can have the form C D XfXxR 
with Xf the coupling constant. The s-channel annihilation cross section is a p-wave 
process which is suppressed by the cross section is given by 


(TUrel 


(s - m|)2 + 2m^ 


( 21 ) 


Similarly, the expression for J/ is 
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Using the (5-function approximation, one finds 


Jf ^ Jf = -r 

^ ^ 8ml 


'X 


(23) 


Again this approximation does not apply to the case with e^j < 0 and |ej^| approaching 
zero. The integration of Eq. (22) dominates in the region near the point = ^R if 
7 ij <C 1, and the integrand decreases rapidly with leaving away from e^. Since the 
situation we considered is near the resonance point (e/j ~ 0), the integral of Eq. (14) 
can be done in the region 0 < yT7nrei/2 < 1. Using the Taylor expansion of the 
complementary error function 


erfc(yT7t>rei/2) = 1 
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and retaining the first order term of Wrei in the series, Jj can be approximated with a 
reasonable precision as follows. 
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b — \/a + e/{, 


c — \/a — e/j, 


d = 




(26) 



Figure 1; Analytical approximations J/, J/ and numerical result JJ? near the resonance 
point when x is a scalar (a) or fermionic (b) DM particle in a specific case. The parameter 
values are taken as /i = 1 GeV, A/ = 1, mij=200 GeV, y^=10, Xf = 20, and F/j = 
0.001 GeV. 


In Figure 1(b), we show the differences in approximate analytical results J/, J/, and 
the numerical result J/ in a specihc case where the parameters are taken as A/ = 1, 
mR=200 GeV, = 10, x/ = 20, and Fjj = 0.001 GeV. As can be seen from the figure. 




































if < 100 GeV, the analytical result agrees with the numerical result JJ well, the 
relative error is less than 5%, and the analytical Jj can obtain the same precision if 
3> 7 i? is satisfied. If > 100 GeV, is no longer valid, but still agrees with 
near the resonance point well and the relative error is within 11 % (in the region 
100 GeV < < 100.01 GeV). On the resonance point, the relic density is given by 


^ 2.755 X 10^ 


45 


8ml 

X. 


'^9* XjMpan ln(2x^^/%^^) - 4/0? 


^GeV 


-1 


(27) 


Unlike the s-wave case, the relic density is inversely proportional to 7 r. 

If X is Majorana fermion, we can write down the Lagrangian of y interacting with 
7? as: £ D the expression for aUrei and relic density is identical to the case of 

Dirac DM. 


3 DM direct detection 


Direct detection experiments search for the signal of DM via their interactions with 
nucleus (for a review, see e.g. [59]). A DM particle can interact with nuclei through 
t-channel scalar R exchange. Since R mixes with the Higgs boson, it can couple to the 
SM fermions with coupling constant rrif/vs, where m/ is the fermion mass and Vs is a 
mass scale parameter. The decay width T^j of the scalar R is given by 


rK = 


E 


ncrjrrij {mj^ — Amj) 2 

STTVg m\ 


(28) 


where ric = 3(1) is the number of color for quarks (leptons), 7 = 1(1/2) for (in)distinguishable 
hnal particles. If y is a scalar DM particle, the spin-independent DM-nucleus elastic 
scattering cross section is given by [60] 


a 


SI 

n 


47rm^m|j 


J ni 


(29) 


where is the DM-nucleon reduced mass + rrin) with is the 

target nucleus mass. /„ stands for the coupling between R and nucleus, which is given 
by 


fn _ V ^ f(n) ^ p{n) \ ^ Rq 

~ ^ JTq r -^JtG / . ’ 

^ rrir, 27 rrir, 


rria 


(30) 


where = 0.011, = 0.0273, and = 0.0447 [61]. The coupling between 

DM and gluons from heavy quark loops is obtained from = 1 — '^q=u,d,sf^q ■, which 
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( 31 ) 


leads to f!j^Q = 0.917. In this case = rriq/vs, then 


gj 0.02056/i^„m„ ^ 


Making nse of Eq. (27), (28), and the latest experimental observation VLji^ = 0.1199 ± 
0.0027 [62], we obtain the expression of for the DM annihilation into SM fermions 
throngh the resonant state R 


4.12 X 10 

^ Y, ncrjrnj{ml-mj)^ 


(32) 


The above expression shows that is proportional to 7 /?. 

For the fermionic DM particle described in section 2.2, it also interacts with nn- 
clei throngh t-channel scalar R exchange. The spin-independent DM-nencleus elastic 
scattering cross section is 


^gj _ 0.08224p2^m2 


o-„ = 


TTVtm 


R 


1.65 X 10 


In(2a: 






m. 


(33) 


Compared with the scalar DM case, is not sensitive to 7 ij. 

For DM direct detection experiments, there is an irredncible backgronnd created 
by the coherent scattering of cosmic neutrinos off target neuclei. This irreducible back¬ 
ground is very difficult to be distinguished from the interactions between DM-nuclei scat¬ 
tering, and it can set a limit on the sensitivity of DM direct detection experiments. Due to 
the neutrino background, the sensitivity of the spin-independent DM-nucleus scattering 
cross section of DM direct detection experiments is limited to 10“"^® cm^ ~ 10“"^® cm^, 
depending on the DM mass [63,64]. If is below the neutrino background, the signal 
of DM particles can not be reached by DM direct detection experiments. 

In Figure 2, we show the relation between and m^ on the resonance point with 
different values of when x is a scalar (a) or fermion (b). In Figure 2(a), for a typical 
DM particle mass ~ 50 GeV, is above the neutrino background when the condition 
7 ij > 2.2 X 10“'^ is satished. In Figure 2(b), is always above the neutrino background, 
as long as the DM particle is lighter than the top quark. 
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Figure 2: The relation between and with different for % being a scalar (a) or 
fermion (b) DM particle. We take the parameters value = 10, x/ = 20. 


4 Resonant annihilation in the real singlet dark mat¬ 
ter model 


In this section, we consider the resonant annihilation of DM particle in the real singlet 
DM model [17-26]. The Lagrangian of the real singlet DM model is [17,19] 

C = CsM + \d,Dd^^D - - XhD^H^H, (34) 

where is the Lagrangian of SM, H is the SM Higgs doublet. The linear and cubic 
terms are forbidden due to a discrete Z 2 symmetry D —)■ —D. D has a vanishing 
vacuum expectation value (VEV) to ensure the DM stability. describes the DM self¬ 
integration strength which is independent of the DM annihilation. It is clear that the 
DM-Higgs coupling Xh is the only one free parameter to regulate the DM annihilation. 
After the spontaneous symmetry breaking, one can obtain the DM mass m|, = rriQ + XHVQ 
with the vacuum expectation value Vq = 246 GeV. In the real singlet dark matter model, 
the DM annihilation cross section is given by 


(JUj-el 


8A>g 

(s — mjjy + 2mD 


(35) 


where T^^ is the total decay width of Higgs which may decay to fermion pairs, gauge 
boson pairs and the real singlet DM pairs if rriH > 2m d [26], the value of r//(-\/s) is 
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given by 
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A / _ 4m^ 12m^ \ - 4m|, 

IGttoq V g V / Svr s 


( 36 ) 


and F/f = r//(\/i) Eq. (35) can be written as 


(TWi-el 


+ 'll' 


where 


(^H 


rri^j — 4m|) 

ml 


and '-)H 


mtrg(mg) 

ml 


(37) 


(38) 


The real singlet DM model is a specihc example for the case of s-wave annihilation 
we have discussed, the relic density for the real singlet DM and near the resonance 
point is analogous with Eq. (19) and Eq. (31), and they can be obtained by substituting 
the parameters ju, jn, and by 2XhVq, ■jh, and en- Figure 3 shows the numerical 
and analytical value of cx®^, and the upper limits for the spin-independent DM-nucleus 
cross section from LUX [65] and XENONIOO [66]. In the hgure, we hnd is above the 
neutrino background and it is not excluded by the result of LUX and XENONIOO near 
the resonance point {mn = mu 12). 

Currently the strongest upper limits on are given by LUX experiment [65] and 
the next generation of DM direct detection experiments can push the upper bound on 
(X®^ down to ~ 10“^’^ cm^ [67]. As the direct detection experiments at present can not 
measure the below 10“^® cm^, so we are uncertain of the existing of the real singlet 
DM near the resonance point. If the future direct detection experiments prove the region 
near the resonance is excluded, the real singlet DM will be removed from dark matter 
candidates. In conclusion, if we want to test the singlet dark matter model thoroughly 
by direct detection, the experiments’ ability should reach the minimum value of 
which is about 2.6 x 10“^® cm^. 


5 Conclusion 

In summary, we have presented an approximate analytical expression for the DM 
relic density of a scalar or fermionic DM particle which annihilates through an s-channel 
scalar resonance which has a narrow decay width. Based on the expression, we have 
investigated the condition for the DM-nuclei scattering cross section to be above the 
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Figure 3: The value of near the resonance point. Parameter values are x/ = 20 and 
= 10 . 

neutrino background. It is found that in Higgs-portal type models, for DM particles 
with s-wave annihilation, the spin-independent DM-nucleus scattering cross section is 
proportional to For a typical DM particle mass ~ 50 GeV, the condition leads 
fo 7i? > 0(10“"^). In p-wave annihilation case, the spin-independent scattering cross 
section is insensitive to and is always above the neutrino background, as long as the 
DM particle is lighter than the top quark. In the real singlet DM model, is always 
above the neutrino background. In order to cover the full parameter space of this model, 
the required sensitivity should reach ~ 2.6 x 10“^® cm^ for the next generation direct 
detection experiments. 
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